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The radiation pressure coupling with vacuum fluctuations gives rise to energy damping and 
decoherence of an oscillating particle. Both effects result from the emission of pairs of photons, 
a quantum effect related to the fluctuations of the Casimir force. We discuss different alternative 
methods for the computation of the decoherence time scale. We take the example of a spherical 
perfectly-reflecting particle, and consider the zero and high temperature limits. We also present 
short general reviews on decoherence and dynamical Casimir effect. 

7— I ■ 

O I. INTRODUCTION AND BRIEF SUMMARY OF DECOHERENCE THEORY 

o : 

The understanding of the quantum-to-classical transition has been the subject of extensive research The core 

of the problem is that the Hilbert space of allowed states of a quantum system is huge, whereas the set of states with 
associated classical properties is a tiny subset of the whole Hilbert space. Some questions that naturally arise are 
the following: which mechanism is responsible for the classical appearance of macroscopic and mesoscopic quantum 
systems? How are those few classical states selected from the huge Hilbert space? The common wisdom is that 

^ — I ■ classicality is an emergent property induced on subsystems by their environment. 

The interaction between a system S and its environment E creates entanglement (i.e. non reducible correlations) 
between the states of the system and those of the environment. Imagine that at a given time (say t = 0) the state 
\^(t)) of S+E is a product state, that is, there are no initial correlations. We have \^>(t = 0)) = |s)|e), where the 
first ket corresponds to the system state (assumed for simplicity to be in a pure state), and the second one to the 
environmental state, also assumed pure. When the two parts begin to interact, entanglement is generally produced. 

CO ' This means that at a later time t, the state | will be given by a linear superposition of the form 

o : 

|*(t)) = | Sl )| £l ) + | S2 )|e 2 >+... (1) 

where {|si)} and {|ej)} are states of the system and environment Hilbert spaces, respectively. If the interaction is 
O ,' such that the states | e^) become approximately orthogonal ((e„|e m ) ~ 5 nm ), then interference between the system 
states \s„) and \s m ) will not be observed. These set of states {|si)} usually have classical properties. Any quantum 
superposition of them is a non-classical state, and quickly decays away into a statistical mixture of the states. The 
coherence of the phase relation between the components of the superposition is lost, and this process is accordingly 
^ known as decoherence. In other words, the environment monitors the different classical alternatives for the system (the 
different states |si)), thereby providing which-way information, even though such information is usually unaccessible 
to the observer. The set of states {|si)} are called pointer states ||, and they are the states within the huge Hilbert 
space of the system that become less entangled with the environment. Perfect pointer states are those that produce 
no entanglement at all, so that an initial product state of S+E will remain a product state throughout the interaction 
time, which means that those states are robust and stay unperturbed by the interaction. All this will be illustrated 
in Section 3 in the particular case where the environment is the radiation field at zero temperature (vacuum field) 
and radiation pressure is responsible for the coupling between system (a mirror) and the environment. 

A possible method to identify pointer states is called the 'predictability sieve criterion' which is based on the 
fact that pointer states are the ones that produce least entropy and remain most pure. Let us explain these concepts. 
The evolution of the closed combined system S+E is unitary, so that the purity of the whole state \^) is preserved, 
i.e. P(ps+e) = Ttp's+e = 1 f° r au times. However, the purity of S is not preserved. To show it one needs to calculate 
the reduced density matrix of the subsystem S by tracing out the environmental degrees of freedom, ps = Tr ePs+e, 
and then P(ps) < 1. The loss of purity can be associated with a loss of information about the system state. When no 
measurement involving the environment is made, the density matrix ps contains the state of knowledge of an observer 
about the system, and purity is a measure of that knowledge. Initially, there is full knowledge of the system state, 
which is described by a single ket state. Subsequently, the interaction with the environment produces entanglement, 
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and part of the information about the system S is lost to the environment, causing a decrease in the purity of the 
system. One can also measure the information content of ps through the von Neumann entropy, S(ps) = — Trps log ps- 
Initially one has full information, and entropy is identically zero; as time goes on, information is lost and entropy 
is produced. Pointer states arc least affected by the environment, so their information content is preserved and 
hence they produce the least entropy. The idea is then to take every state in the Hilbert space of S, calculate the 
von Neumann entropy at time t produced via interaction with the environment starting from the given state of the 
system, and order the initial states in a tower of increasing entropy. The states that lie at the bottom of such a tower 
are candidates for pointer states. Finally, one must check that those states remain at the bottom when the time t 
when entropy is calculated is changed, so that the states are robust pointers. If one can satisfy these conditions, one 
gets the pointer states of the system S. If not, there are no pointer states for the problem. When these ideas are put 
into practice, it is more convenient to work with the linear entropy, defined as S{ps) = 1 — P{ps) = 1 — Trp|. 

The determination of pointer states for a given problem depends both on the system and environment dynamics, as 
well as on their interaction. There are three different regimes: 1) When the system Hamiltonian is irrelevant, pointer 
states are given by the eigenstates of the system's operators contained in the interaction Hamiltonian. A typical 
example is that of a measuring apparatus (the system) that has no internal dynamics, measuring an external reservoir 
(the environment), say a photocurrent. 2) When the system dynamics as well as the interaction are relevant, pointers 
come from an interplay between the two. The most thoroughly studied example is that of quantum Brownian motion 
(QBM), in which a particle is coupled to a set of harmonic oscillators in a thermal state, thereby suffering decoherence 
p] . Although the coupling between the system and the environment is of the type position-position, pointer states 
are not position eigenstates of the system because the self dynamics of the harmonic oscillator interchanges position 
and momentum every quarter of a cycle. It turns out that the interplay between interaction and self dynamics leads 
to coherent states as pointers ||. This second case is also the relevant one for this paper. 3) Finally, when the 
environment evolves much slower than the system, pointers may correspond to energy eigenstates of the system's 
Hamiltonian H. 

In the above we have ignored the information contained in the environmental state, and that is the reason why one 
traces over the environmental degrees of freedom in order to find the reduced density matrix of the system. However, 
the information lost to the environment could be, in principle, intercepted and recovered. Performing measurements 
on the environment one may extract information about the decohering system. In || it is shown that the preferred 
pointer states remain unchanged, even when that information is kept and modifies the dynamical evolution of the 
system. 

Another related way to study the dynamical process underlying in the quantum-to-classical transition is via phase 
space representations of the reduced quantum dynamics for the system S. Among the many possible representations, 
one of particular interest is the Wigner function W(x,p), which is defined as a Fourier transform of the reduced density 
matrix. W(x,p) is a pseudo probability distribution in phase space, and encapsulates the quantum coherence of the 
system in interference fringes that take both positive and negative values. Imagine that one starts with an initial state 
for S which is highly non classical, such as a cat state |cat) = l/y2(\a) + | — a)), where \a) is a coherent state with 
large amplitude (|a| 3> 1). The corresponding Wigner function will have interference fringes, showing the quantum 
nature of the state. However, when the system is put in contact with the environment and each component of the 
state becomes entangled with almost orthogonal states of the environment, the interference fringes will be washed 
out. In the end the Wigner function becomes positive defined, with two peaks at the values corresponding to ±a, 
as a true probability distribution. Decoherence transforms the initial pure state into a mixture of the two coherent 
states | ± a) . 

Until not very long ago the ideas of decoherence were restricted to the theoretical domain. Recent experimental 
developments have succeeded in studying in real time the process of decoherence in the laboratory, and have tested the 
predictions of the theory. Here we shall mention a few experiments that have been a hallmark in those developments. 
First, in the field of cavity QED, superposition states of photons and Rydberg atoms have been created within 
high-Q microwave cavities. Cat states of around 3 photons have also been produced, and it has been studied how 
they decay due to decoherence [||. The coherence of the state was monitored with the help of a measurement of 
correlations between two consecutive atoms crossing the cavity EjJ. Second, in the field of ion trapping, methods for 
creating superposed motional states of ions were developed, as well as schemes of environment engineering to protect 
those states from decoherence fLl[ . Finally, it has been possible to push the size of the cat states further into the 
macroscopic realm by generating a mesoscopic cat inside a rf-SQUID. The two components of the cat correspond to 
superconducting currents moving either clockwise or counterclockwise, each containing around 10 9 Cooper pairs |L2] . 

In what we have discussed so far decoherence has been portraited as a "good" effect, in the sense that is responsible 
for the quantum-classical transition and the appearance of our classical world. Decoherence can also have a "bad" 
role in the field of quantum computation and quantum information processing. There one performs logical operations 
making use of the superposition states of quantum mechanics. For such operations to be successful it is very important 
to maintain the relative phase between the components of the superpositions all along the operations. If decoherence 
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acts, it produces quantum errors that must be somehow corrected. Several methods have been proposed to minimize 
the effects of decoherence (see |0 for an example in nuclear magnetic resonance) . 

The prototype calculation of environment induced decoherence is the heuristic position-position interaction Hamil- 
tonian for describing quantum Brownian motion, where the environment is taken to be a collection of harmonic 
oscillators. Although such a model is quite useful for studying many physical processes associated to dissipation 
and decoherence of a quantum system, the results that follow from it do not apply to every situation. That is, it 
is necessary to perform a case by case analysis in order to compute physical observables, such as decoherence and 
damping rates, how they scale with the parameters of the system, the environment, and their coupling, etc. For usual 
environments (thermal atoms, thermal light, phonons, etc.) it is in principle possible to design engineering schemes 
to protect the state of the system from decoherence, for example by reducing the coupling to the environment. 

Then, the following question naturally arises: is it possible, at least in principle, to have arbitrarily weak deco- 
herence? In this paper, we consider a fundamental source of decoherence that cannot be 'turned off': the radiation 
pressure coupling with the vacuum field pi] ]. As reviewed in Sec. 2, any particle not completely transparent unavoid- 
ably scatters vacuum field fluctuations. This type of coupling is responsible for the Casimir effect. More generally, 
photons are created out of the vacuum field when moving boundaries are considered, an effect known as dynamical 
Casimir effect or motion-induced radiation. In Sec. 3, we show how the dynamical Casimir effect engenders decoher- 
ence. Our emphasis is on the basic physical ideas, and most of the calculations are referred to |l5f| , but we also briefly 
discuss a model alternative to the one employed in this reference. 



II. DYNAMICAL CASIMIR EFFECT 

The Casimir effect is perhaps the simplest and most striking effect of the quantum vacuum field (see for 
reviews). The essential idea is that the boundary conditions modify the spectrum of the radiation field, and thereby 
its zero-point energy. This modification has direct physical consequences, leading, for example, to an attractive force 
between two parallel perfectly-reflecting plates (of surfaces A) and at a distance L, given by 

240 L 4 

A series of recent experiments Jl9| reported precision measurements of the Casimir force in agreement with the predic- 
tions of Quantum Electrodynamics, although more complete theoretical calculations, taking into account corrections 
due to finite temperature and conductivity as well as to roughness and geometry of the surfaces are partially yet to 
be done @. 

The Casimir force may also be computed by taking the average of the Maxwell stress tensor over the field vacuum 
state H^]. This method suggests that the Casimir force is itself a fluctuating quantity, as noted by Barton. Its 
fluctuations were first computed for plane perfectly reflecting mirrors [^2), and later for spherical and spheroidal 
particles [^3| . More generally, any particle scattering the radiation field is under the action of a fluctuating radiation 
pressure force exerted by the vacuum field, even in the situations where the average force vanishes (for example a single 
plane mirror at rest). The coupling responsible for those fluctuations also gives rise to a dissipative force, when the 
particle is moving in vacuum. Dissipation of the mirror's mechanical energy is needed to enforce energy conservation, 
since the motion induces the emission of pairs of photons (for reviews see p7| ]2§|]). Because of their common physical 
origin, fluctuations and dissipation are related by a very general result ]29[, whose most known application is the 
Einstein relation between diffusion and friction coefficients for a Brownian particle in the high-temperature limit. 
This connection provides a very useful tool for deriving the response to an external small perturbation from the 
fluctuations in the unperturbed case. Linear response theory was employed by Jaekel and Reynaud to infer the 
vacuum radiation pressure force on partially-reflecting moving mirrors in the one-dimensional (ID) case fjCfl . For a 
single perfect mirror (position x(t)) the force is given by 

6nc 2 dt 3 ' 

a result first obtained by solving the boundary conditions of a moving mirror in the long wavelength approximation, 
and assuming the effect of the motion to be a small perturbation |3l[] . Eq. (Q) was also derived as the n — > oo limit 
of a moving half-space of refractive index n f32|| . It also corresponds to the nonrelativistic approximation of the 
exact result (for a perfect mirror) derived with the help of a conformal coordinate transformation to the co-moving 
frame @. 

Since the wave equation in three dimensions is not invariant under a general conformal transformation, only ap- 
proximated methods are used in this case. The dissipative force on a plane mirror was computed within the long 
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wavelength approximation for a scalar and electromagnetic |34| field models. The an gula r and frequency dis- 



tributions of the emitted radiation were also computed for a single plane moving mirror 35 1, a moving dielectric 



half-space |36| |37| and two parallel plane mirrors |38|]. Linear response theory was employed to derive the dissipative 
force on moving spheres |3j| . Small but otherwise arbitrary time-dependent deformations of an initially plane surface 
were analyzed with the help of different approaches: linear response theory E3|, long wavelength approximation fl4l| ], 
and path integrals Q . 

The magnitude of the dynamical Casimir effect may be illustrated with the following example, which we shall 
discuss in detail in Sec. 3. We consider that the 'mirror' is a particle of mass M in a ID harmonic potential, such 
that the oscillation frequency is ujq. From Eq. (Q), the equation of motion reads 

d x 2 h d x , . 

For any situation of physical interest, the zero point energy is much smaller than the rest mass energy: huio <C Mc 2 . 
In this case, (||) has solutions corresponding to oscillations damped at the rate 

r = ^V, « W0 , (4) 



UirMc 2 



showing that the dynamical Casimir effect provides a tiny perturbation of the free oscillations. 

As could be expected, a larger effect takes place when field modes of a cavity resonator are coupled to the moving 
boundaries, mainly when the mechanical frequency lies close to a given cavity eigenfrequency. Moore considered a 
scalar ID field inside a cavity where one of the mirrors follows a prescribed motion The field modes were formally 
built in terms of the solution of a functional equation. This method was later developed |Q and extended to the 
case where the two mirrors are set in motion [|45|. The case of partially-transmitting mirrors was also calculated, 
allowing for a reliable estimation of the orders of magnitude for the rate of transmitted photons and the number of 
photons inside the cavity at steady-state [Q. So far, few three-dimensional (3D) calculations along these lines have 
been reported. A rectangular cavity made of perfectly-reflecting moving mirrors |47[ p8[ , and a spherical bubble with 
time-dependent radius [|49| were analyzed, the latter motivated by the problem of sonoluminescence. 

In this article we only consider a single scatterer, so that no resonant enhancement takes place. In this section, 
we have shown that the radiation pressure coupling gives rise to energy damping of a particle scattering vacuum 
fluctuations. In the next section, we show that it also destroys the quantum coherence of the particle. 



III. DECOHERENCE AND THE CASIMIR EFFECT 



Most treatments of the dynamical Casimir effect consider the particle that scatters the vacuum field (the 'mirror') 
to follow a prescribed motion (an exception is Ref. ]5Q| , which considers fluctuations of position of a particle driven 
by vacuum radiation pressure). In this article, however, we want to focus on the particle as the dynamical degree of 
freedom of interest. More specifically, we analyze how the radiation pressure coupling destroys the quantum coherence 
of an initial superposition state of the particle. 

We consider as before that the particle is in a harmonic potential well, corresponding to a frequency of oscillation 
uio. The connection with the previous approaches, where the (classical) particle is assumed to follow a prescribed 
oscillation, is made by taking a coherent quasi-classical state \a) for the particle, so that the combined particle-field 
state at t = is 



|*(f = 0)> = |a>|0>, (5) 

where we have assumed that the field is initially in the vacuum state |0). The oscillation gives rise to the emission of 
photon pairs at time t at the field modes Ai and A2, with probability amplitudes b{\\, A2, t) : 

|*(i)) = \a) ( B(t)\0) + J2 b(X 1 ,X 2 ,t)\X 1 , A 2 ) I , (6) 

\ Ai,A 2 / 

where B(t) is such that this state is normalized. As discussed in Sec. 2 (see, in particular, Eq. (||)), the energy damping 
associated to the dynamical Casimir effect is very small. This effect, and more generally the recoil of the particle, is 
neglected in (^|), where the particle state is assumed not to be modified. Even at this level of approximation, there 
is decoherence, as we show by taking the initial state of the particle to be the cat state |cat) = (|a) + | — a))/^/2, 
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an example already mentioned in Sec. 1. It corresponds to the coherent superposition of two wavepackets oscillating 
out-of-phase in the harmonic potential well. The amplitudes 6(Ai,A2,£) depend on the phase of the oscillation, so 
that they have an opposite sign when we take the state \ — a). Since the evolution operator is linear, the complete 
state at time t is the superposition of the r.-h.-s. of (||) with the analogous state for | — a). It turns out to be an 
entangled state of the form discussed in (|l]): 

\*(t)) = \a)\eM(t)) + \-a)\e^{t)), (7) 

with |e^(t)) = B(t)\0) ± J2\ ± \ 2 &( Ai , Aa , *) | Ai , A2). These field states work as tags for the particle states, providing 
which-way information about the phase of the oscillation. As time goes on, the information gets better defined, since 

(e(->(t)|e< + >(t)) = l-2£ \b(X x ,X 2 ,t)\ 2 (8) 

Ai,A2 

decreases as the probability for photon emission increases. When the emitted photons are not detected, all the relevant 
information about the particle is contained in the reduced matrix p(t) = Trp (| ^(i)) (VP (t) |), where the trace is taken 
over the field states. Since the interference term is gradually washed out as a consequence of the photon emission effect 
and the corresponding entanglement with the field, p(t) decays into the statistical mixture p m = (\a)(a\ + \—a) (— a|)/2. 
The corresponding time scale f <j may be computed [jl5| from Eq. (^) , and turns out to be proportional to the energy 
damping time 1/T, which is related to the two-photon probabilities by energy conservation: 

Eq. (JoJ) also holds when the coupling with the environment is described by a heuristic master equation in the Lindblad 
form (derived with the help of the rotating- wave approximation) |5l[ , as well as in the case of position-position 
coupling to a zero-temperature environment of harmonic oscillators, and has a very simple interpretation fll6| ]: if 1/r 
is the time needed to damp the energy 2|a| 2 ?iwo, it corresponds to the emission of 2\a\ 2 pairs of photons (each pair 
has a total energy equal to hu)a). On the other hand, coherence is much more delicate than energy, since a single 
photon provides which-way information that destroys the quantum phase of the cat state. Hence the decoherence 
time is the time scale for the emission of a single photon. Since 4|a| 2 photons are emitted during the time interval 
1/r, the time for a single photon scales as in the r.-h.-s. of @. 

Eq. only holds when \a\ 3> 1. In this limit, decoherence is much faster than damping, justifying the approach 
of neglecting the decay of the amplitude a of the coherent states in (||) and (^). This is of course in line with the 
idea that in the 'macroscopic' limit weird non-classical states are extremely fragile and difficult to observe. For truly 
macroscopic systems td is so short that no experimental monitoring of the decoherence process is possible. However 
the validity of Eq. @ is restricted by the additional condition that decoherence is slower than the free oscillation 
(this condition is fulfilled by the experiments || discussed in Sec. 1). In this regime, the particle oscillates several 
times in the potential well before coherence is lost, and the r.-h.-s. of Eq. (|^) may be written in terms of the distance 
Aa; = 2^/2Ti/Mwo\a\ between the two wavepackets when they are at their turning points (M is the mass of the 
particle): 

where Aieo = y/K) {2Mloq) is the position uncertainty of the oscillator ground state. Eq. (|Io|) shows more explicitly 
that the decoherence rate scales as the squared distance in phase space between the two components of the cat state. 
In Eq. (^|), the distance is expressed in terms of the squared difference Aa = 2a between the amplitudes of the two 
coherent states | ± a). Such dependence, already experimentally observed in was fully verified in [fi"lf . Thus, the 
decoherence rate is directly connected to the quality of which-way information, for the possibility of resolving the two 
wavepackets is quantified by the distance between them divided by their width Axo- 

The second factor entering in the r.-h.-s. of (||) is the damping coefficient T. Rather than a phenomenological 
constant, here T quantifies the strength of the radiation pressure coupling to the vacuum field, and is calculated from 
first principles. As discussed in Sec. 2, it may be obtained directly from the expression for the dissipative radiation 
pressure force on the particle. In the ID case, T is given by Eq. ([IJ), which jointly with Eq. (^) yields 

3 2?r 

(y/cY ujq 
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where v — y/2fuvo/M \a\ is the velocity of the wavepackets at the moment they cross the bottom of the potential well. 
Therefore, in the nonrelativistic limit considered in this paper, decoherence is much slower than the free oscillation. 
The ratio between the two time scales is even larger when considering the real 3D case. If we take a spherical perfectly- 
reflecting particle of radius R smaller than the range of oscillation, then loqR/c < v/c <C 1. Since the relevant field 
modes have frequencies of the order of u)q, in this limit the particle is much smaller than typical wavelengths (Rayleigh 
scattering regime), and hence is weakly coupled to the field. The dissipative force in this regime was calculated in [p9| ; 
the resulting damping coefficient scales as the squared polarizability of the sphere, leading to an additional factor 
(^ i?/c) 6 : 

324 / c \ 6 2vr /c\8 2tt 

td = 7T« — p — » - — ■ (12) 

It is also possible to analyze the decoherence effect in a more complete theoretical framework, where the dynamical 
radiation pressure coupling between particle and field is fully taken into account. This approach also accounts 
properly for damping of the particle's energy, as well as for additional effects resulting from the coupling with the 
field. Moreover, it allows us to analyze decoherence in the more general case of an arbitrary temperature of the field. 
An ab-initio Hamiltonian model for the particle-field system was derived from first principles in Ref . [^6| . This model 
was the starting point for the discussion of decoherence in Refs. |Q and p5[ . The field scattering corresponds to 
frequency dependent reflection and transmission coefficients that satisfy the passivity requirements discussed in [p2[ . 
This means that the dynamics of the particle does not suffer from the instabilities associated to the model of a perfect 
mirror (as well known from classical electron theory, Eq. (||) is plagued with 'runaway' solutions). 

Here we describe the radiation pressure coupling with the alternative, more intuitive model, where the interaction 
Hamiltonian corresponds to the energy transfer between field and particle: 

H int = -xF, (13) 

where F is the radiation pressure force on the particle, and x its position. This type of model was extensively employed 
in several contexts associated to the dynamical Casimir effect ^7j. Here we focus on the limit where the particle 
perfectly reflects the (ID) field, but a discussion of partially-reflecting mirrors along these lines is also possible. As 
shown below, it leads to results for the decoherence and damping rates in agreement with those found in Ref. fl5|j . 

Starting from (|l3|), we derive a master equation for the reduced density matrix of the particle. It is similar to 
the master equation for QBM, derived from the position-position interaction Hamiltonian. Technically, the essential 
difference arises from the fact that the force operator F in (jl^) is quadratic in the field operators, which leads to a 
damping coefficient that depends on the state (and hence temperature) of the field (reservoir). Although the formalism 
relies on a ID model, the final results may be generalized to the 3D case. 

We write the master equation in terms of the Wigner function W(x,p, t) : 

v o d 2 W d 2 W 

dtW = -Yjd x W + Mu* 2 xd p W + 2Td p {pW) + D 1 ^— - D 2 ——. (14) 
M op z oxop 



The first two terms in (14) correspond to the harmonic oscillation in the potential well, with a frequency ujq = ujq + Slu 
modified by the coupling with the field (on the other hand, when the interaction Hamiltonian is linear in the momentum 
of the particle, a mass correction appears J2(|). The remaining terms describe non-unitary evolution. The damping 
as well as the diffusion coefficients D\ and D 2 are time dependent and given in terms of correlation functions of the 
force operator. The diffusion coefficients are related to the symmetric correlation function: 

a FF (t) = ({F(t),F(0)}), (15) 

where the brackets denote the anticommutator, and the average is taken over the field state (thermal equilibrium, 
temperature T). 

The term proportional to D 2 in ( |l4| ) yields a negligible contribution, so that we focus on D\ : 

1 f duj 

D t(t) = 2 / 2~°>FMsync t M, (16) 
where a^f [w] is the Fourier transform of <JFF{t) and 

sync t (w) = 

LO — U)Q 



G 



is a function peaked around u> = ujq of width 2%/t. Clearly, for a time t long enough, the function sync t (u>) is so sharply 
peaked that <jff[lu] is approximately constant over the short frequency interval that contributes in the integral in 
Eq. (pif), and hence may be replaced by its value at ui — loq. In this case, we find 

Di(t —f oo) = jCTff[wo]. (17) 



A sufficient (and also necessary at T = 0) condition for the validity of (|17|) is u>ot ^ 1. In other words, for times 
much longer than the period of oscillation, the field fluctuations at frequency u>o provide the dominant contribution 
to diffusion. 

The damping coefficient is likewise connected to the average value of the commutator of the force operator taken at 
different times (anti-symmetric correlation function). When the interaction Hamiltonian is linear in the operators of 
the environment, as in the position-position model, the commutator is a c-number times a delta function (in time), and 
as a consequence, the damping coefficient has a constant value that does not depend on the state of the environment. 
As already mentioned, this is not the case for radiation pressure coupling. In particular, the damping coefficient 
depends on the temperature of the field, as could be expected having in mind the Stefan-Boltzmann law. At zero 
temperature, we recover the result given by Eq. (^). 

We calculate the pointer states using the predictability sieve criterion, discarding all information about the envi- 
ronment, as discussed in Sec. 1. We start from the master equation, and evaluate the rate of change of linear entropy, 
assuming an initial pure state. It is straightforward to show that the entropy is minimized for minimum uncertainty 
Gaussian states, hence the pointer states are the coherent states. This result agrees with the well-known fact that 
coherent states provide the closest possible realization of a classical state of oscillation, given the constraint imposed 
by the Heisenberg uncertainty relation. In short, coherent states remain approximately pure because they do not 
entangle with field states, at least for times shorter than the damping time 1/I\ as shown by Eq. (^|). 

On the opposite extreme in Hilbert space, superpositions of coherent states are highly nonclassical and cannot last 
when the distance between the two components is large. This may be analyzed in detail from Eq. (fhf). The coherence 
of the initial state is imprinted on the Wigner function in the form of an interference term W m t that oscillates in 
phase space. When the two state components are spatially separated by a distance Ax, the oscillation is along the 
axis of momentum: Wi n t(a;,p) ~ cos(Axp/h). Thus, according to Eq. (|l4|), diffusion washes out this oscillatory term, 
the faster the larger the value of Ax. With an additional factor of 2 to take into account the average over several free 
rotations of the state in phase space, we find 

^ = 2 a|*F (18) 

To derive the decoherence time when the field is in the vacuum sate, we compute the correlation function o~pf[ujq\ 
at zero temperature. When replacing the result for D\ as given by ( |l7|) into (|l|), we obtain the same result already 
derived in this Section by a more elementary method. 

For finite temperatures, the spectrum is approximately constant at low frequencies, so that ( |l7|) also holds when 
u>o <C kT/h (k is the Boltzmann constant), including the free particle limit ujq = 0, provided that the entire frequency 
interval around loq is contained in the low frequency part of the spectrum, which corresponds to the condition 
1/t -C kT/h. The damping coefficient may calculated in the high temperature limit as well, and the results are in 
agreement with Einstein relation 

D l = 2MkTY. (19) 

More generally, we may derive a relation between diffusion and damping coefficients valid for arbitrary values of tem- 
perature Jig ], including T = 0, starting from the general relation between symmetric and anti-symmetric correlation 
functions (fluctuation-dissipation theorem). 

The decoherence time for high T is derived by replacing ( |l9| ) into (|l^). As we discuss below, usually in this limit 
decoherence is faster than the free oscillation, so that, contrary to the T = case, there is no average over many 
oscillations in this case. To describe the decoherence process, we must evaluate the diffusion coefficient at a time t 



much shorter than td- Hence, we are allowed to use its asymptotic value as given by the Einstein relation (19) only if 



we assume that td ^> h/{kT). The resulting expression is very general B, and also holds in the free particle case: 



where At = hj \J1MkT is the de Broglie wavelength of a particle of mass M in thermal equilibrium. Eq. (EG) has 



a form similar to (10), except that now the reference of distance is set by thermal fluctuations instead of zero point 
fluctuations. 
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In order to complete the evaluation of the decoherence time, we need to evaluate the damping coefficient T in the 
high temperature limit. We consider as before a sphere of radius i?, which is usually much larger than typical field 
wavelengths, which are of the order of hc/(kT) (except for very low temperatures or very small spheres). In this 
short-wavelength regime, the radiation pressure force may be calculated by replacing the surface of the sphere by a 
collection of tangent planes, and the final result reads 

47T 3 (fcT) 4 2 dx 

The force scales with the surface or cross section of the sphere, and is proportional to T 4 , in agreement with Stcfan- 
Boltzmann law. As opposed to the vacuum case, here we have a true friction force, i.e. proportional to the velocity 
of the particle and not to higher-order time derivative as in Eq. (||) (the thermal field is not Lorentz invariant). 

In the free case (u>q = 0), T is simply the coefficient multiplying the velocity in Eq. (^l]) divided by M. Then, with 
the help of @ we find 

_ 45 h 5 c 4 

td ~ 8tt 3 {kTfR 2 (Axf ( ' 

Eq. (B2) shows that the decoherence time depends strongly on temperature (the same temperature dependence was 
found in Ref. pi)] ). Even at the temperature corresponding to the cosmic background radiation, T = 2.7K, radiation 
pressure is a very efficient source of decoherence. As an example, for R = 1cm, we have ta = 2.7 x lCP 21 /(Ax[m]) 2 s, 
which is in the nanosecond range for a separation Ax = lyum. 



IV. CONCLUSION 



The master equation provides a complete description of the particle dynamics when no measurement on the field 
is made. It accounts for the renormalization of the oscillation frequency, damping, and diffusion and the associated 
decoherence effect. It also allows for the determination of the pointer states, and all that for any temperature T. On 
the other hand, the decoherence time scale at T = may be calculated by a simpler approach, in which wc follow the 
evolution of the complete particle-field state, calculated with the help of the superposition principle, and trace over 
the field at the very end. This approach explicitly shows that decoherence results from entanglement between particle 
and field states. 

The decoherence induced by radiation pressure coupling with vacuum fluctuations is a very slow effect, when 
compared with the the time scale of the free evolution. Yet, it is remarkable, from a conceptual point-of-view, that 
classical behavior of a macroscopic system emerges from the formalism of Quantum Mechanics itself, even though in 
very long time scale, provided that the quantum vacuum radiation field is taken into account. 

P.A.M.N. thanks M.T. Jaekel, A. Lambrecht and S. Reynaud for discussions, and CNPq, PRONEX and FAPERJ 
for partial financial support. 
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